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PBBFAnW 

A  portion  of  th«  earth's  atmosphere  that  has  received  relatively 
little  Investigation  Is  the  mesosphere,  lying  between  20  and  80  km. 
Knowledge  about  a  variety  of  electrical  properties  and  phenomena  of 
this  region  of  the  atmosphere  Is  especially  scanty.  The  mesosphere 
Is  strongly  responsive  to  a  variety  of  excitation  effects.  Including 
bombardment  by  certain  classes  of  solar  protons  and  by  solar  ultra¬ 
violet  and  x-radlatlons .  These  effects  can  produce  notable  variation 
In  the  electrical  parameters  of  the  mesosphere. 

A  previous  RAND  Nemorandum,  RM-3230-PR,  Effects  of  Cosmic 
Radiation  on  the  Extremely  Low  Frequency  Properties  of  the  Mesosphere, 
examined  some  electrical  properties  of  the  mesosphere  that  are  Import¬ 
ant  to  electromagnetic  phenomena  on  the  earth  —  micropulsations, 
extra-low-frequency  noise,  and  hydromagnetlc  waves.  This  Memorandum 
Investigates  the  electromagnetic  structure  of  the  mesosphere  with 
such  geo-electromagnetlc  phenoMna  In  mind.  The  analysis  Indicates 
possible  sources  for  certain  stratification  and  ducting  features  of 
the  mesosphere  that  become  prominent  for  slowly  varying  electrical 


effects 


ABSIRACT 


Electrical  disslpatloa  In  the  mesosphere  is  dominated  by  gas 
particle  collisions.  Hence,  the  electrical  properties  of  this  region 
are  controlled  by  the  density  of  the  atmosphere.  For  low-frequency 
propagation,  the  electromagnetic  properties  of  the  mesosphere  can 
be  described  in  terms  of  variations  with  altitude  of  the  local  complex 
refractive  index.  For  the  simple  condition  of  excitation  of  the 
mesosphere  by  a  single  kind  of  ionizing  radiation,  the  complex 
refractive  index  should  have  exponential-scale  variation  with  altitude 
of  less  than  12  km  for  the  lower  mesosphere.  Sudden  onset  of  more 
complex  conditions  of  excitation  may  ''^Iggef’^^e  mesosphere  into  a 
state  of  '^helvlng^or  ducted  structure  in  its  electrical  properties. 
Depending  on  the  details  of  the  excitation,  such  ducting  will  contain 
regions  of  negative  gradients  in  the  altitude  variation  of  the  electrical 
properties.  Instead  of  merely  the  positive  gradient  with  a  l2''km  upper 


limit. 
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nie  electrical  properties  of  the  lower  Ionosphere  and  of  the 
mesosphere  Influence  several  geo-electromagnetlc  phenomena  In  the 
1-10,000  c/s  portion  of  the  spectrum.  Different  atmospheric  models 
show  a  variety  of  such  phenomena:  micropulsations  (Bomke.  1962), 
cavity  oscillations  In  the  earth— Ionosphere  air  space  fBalser  and 
Wagner.  1962),  and  distant  propagation  of  lightning  fields  (Holzer 
and  Deal.  1956;  and  Holzer.  1958).  The  models  used  Include  sharply 
bounded  Ionospheres  (Schumann,  1952,  1957),  special  distributions  such 
as  Epstein  liters  (Epstein.  1930;  and  Schelkunoff .  1951),  and  two-layer 
models  (Walt.  1960a) .  Exponential  models  of  the  variation  of  the 
electrical  properties  of  the  Ionosphere  with  altitude  have  been  used 
by  Walt  (1960b) . 

In  a  previous  paper  (Goldberg.  1962),  the  author  discussed  the 
extension  of  such  modeling  to  three  layers  of  the  mesosphere:  upper 
(60  km  to  80  km),  middle  (40  km  to  60  km),  and  lower  (20  km  to  40  km). 
The  present  paper  furthers  that  discussion.  It  examines  distinctive 
electromagnetic  characteristics  of  the  mesosphere,  e.g.,  the  conditions 
for  negative  altitude  gradients  of  refractivlty,  those  for  relatively 
large  gradients  of  refractivlty,  and  those  for  abrupt  transitions  or 
"trigger  situations." 


We  will  use  the  complex  Index  of  refraction  as  a  primary  electrical 


parameter  for  describing  the  niesosphere ' s  electrical  structure.  In 
terms  of  the  electric  number  density  n^,  and  the  electron  collision 
frequency  v,  the  complex  Index  Is 


M 
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l+(v/u))^ 
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(U)p/U))^(v/U)) 
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where 

M  >  the  complex  Index  of  refraction, 

(Up  ■  the  angular  plasma  frequency,  (n^e  /m  e^)  , 
e  *  the  electronic  charge, 
m  ■  the  electronic  mass, 

(u  ■  the  angular  frequency  of  the  excitation  field, 
n^  *  the  electron-number  density, 

~  the  absolute  permittivity  of  free  space,  and 
j  -  . 

The  refractive  Index  Is  convenient  for  denoting  electrical  properties 
of  the  mesosphere  In  dimensionless  values.  (In  terras  of  local  values 
of  the  basic  electrical  properties  of  conductivity  a  and  relative 
dielectric  constant  the  complex  Index  Is  given  by  the  relationship 
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III  ■  KHilZATIOH  PRODOCTKat  MECHANISMS  AHD  MODELS 

Over  most  of  the  mesosphere,  the  electron  density's  verletlon 
with  altitude  can  be  represented  by  single  exponentials  or  by 
coaq>oslte8  of  exponentials.  Hie  distribution  of  electron  density  with 
altitude  depends  (among  other  things)  on  the  number  and  type  of  Ionizing 
radiations  present  and  on  various  processes  of  electron  attachment, 
detachment,  and  recombination. 

The  exponential  form  of  representation  Is  developed  here  for  very 
idealized  conditions,  and  methods  for  treating  the  more  complex  cases 
are  discussed  briefly.  The  fact  that  even  the  more  complex  situations 
can  be  represented  fairly  simply  by  exponentials  Is  a  definite  advantage. 

For  the  simplest  condition  ~  a  minimum  of  participating  processes 
In  equilibrium  —  the  local  electron  density  Is  given  by  the  expression 

■ 

where  n^  Is  the  electron  density,  q  Is  the  rate  of  production  of 
electrons  due  to  action  of  a  single  Ionizing  source,  and  the  equality 
of  3  and  ^^(02)  characterizes  the  loss  of  electrons  due  to  attachment, 
with  no  other  mechanism  (such  as  detachment)  being  Important.  When 
mechanisms  like  detachment  occur,  the  electron  density  distribution 
for  much  of  the  mesosphere  can  still  be  represented  by  single  or 
combined  exponential  expressions. 

As  shown  elsewhere  (Goldberg.  1962),  even  a  single,  relatively 
Ineffective  source  of  Ionization  (such  as  cosmic  radiation  at  night) 
produces  refractlvlty  values  In  the  mesosphere  that  are  significant  for 
an  elactromagnetlc  excitation  with  slow  time  variations.  In  the  daytime 
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and  during  solar  storms,  the  electrical  structure  of  the  masos^are  is 
onich  more  complicated,  owing  to  the  action  of  additional  sources  of 
ionisation.  Ihese  more  complex  conditions,  irtiich  are  excluded  from 
this  analysis,  can  be  treated  by  combining  similar  models  over  different 
ranges  of  altitude. 

For  example,  the  electron  density  can  be  represented  by 

n 

e 

where 

Q  ■  the  total  rate  of  production  of  electrons  ***  ^2  ^  ^3  ***  *  *  * 
by  primary  ionising  mechanisms,  1,  2,  3,  etc., 

\  ■  the  ratio  of  negative  ion  to  electron  densities  vT/n^, 

■  the  rate  coefficient  for  dissociative  recombination,  and 
•  the  rate  coefficient  for  Ihompson  three>body  ionic  re¬ 
combination. 

To  compare  this  expression  for  n^  with  the  expression  for  the 
sljiq>lest  case  (n  ■  q/9),  the  above  expression  can  be  put  in  the  form 
of 

where  N  .  is  the  value  of  M  in  the  lisiiting  case  of  H  ■  n  ■  q/0. 

This  allows  N  to  be  expressed  in  the  form  of  a  product  of  three 
exponentials.  The  values  of  H  for  both  the  lower  and  upper  mssosphara 
can  be  represented  quite  well  by  such  a  confined  exponential,  as 
indicated  by  comparison  with  the  values  of  Holer  (1960), 

(1959) ,  Reid  (1961) ,  and  Landmark  and  Lied  (1961) ,  for  various  conditions 


(l-t\)(a^-w^\) 
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of  the  meeoephere.  Ihe  followlog  enalyeis  of  the  eltltude  verletione 
will  be  confined  principally  to  the  case  of  although  It  will 

make  some  use  of  composite  exponential  representations. 

For  the  case  of  n  "  q/3,  an  exponential  representation  of  the 
altitude  variation  of  electron  density  n^  due  to  an  Ionisation  rate  q 
can  be  obtained  by  writing 

“  ‘»oq  » 

idiere  z  Is  the  reference  altitude  at  which  q  ■  q  .  and  where  Z.  is 
oq  oq  q 

the  altitude  scale-factor  for  q.  Similarly,  P  can  be  written  as 

P  “  Pop  “‘p[’f*‘'oP>*p]  » 

where  P^  is  the  value  of  P  at  the  reference  altitude  z^  and  Z^  is 
the  altitude  scale-factor  for  the  exponential  variation  of  p. 

These  expressions  for  q  and  P  give  the  general  form  for  the 
exponential  variation  In  altitude  of  the  electron  density: 

-if)  • 

*“•  “o,?  ■  W'oe  • 

Also,  In  terms  of  composite  variables  to  be  used  later,  the  electron 
density  Is 

“e  "  “eo  “P  <‘on/2pq>  *  "  “eo  “P  *  “P  ' 

where  n  is  the  value  of  the  electron  density  at  a  common  refarance 
eo 

altitude,  z  *  s  a  *  *  ^«te  Z..  Is  the  composite  altitude  scale- 
'  on  op  oq  pq 

factor. 
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and  where  Z  ■  Z.  %ihen  z  ■  z  .  *  > 

n  on  op  oq 

Note  that  the  electron-denalty  altitude-variation  factor,  Z^^, 

can  be  either  positive  or  negative  «dien  the  component  scale-factors, 

Zq  and  Z  ,  are  positive.  This  fact  (and  the  possibility  under  various 
P  Q 

conditions  In  the  mesosidiere  that  the  altitude  scale-factors  Z^  and  Z^ 
can  vary  considerably  In  their  magnitudes  as  well  as  In  their  signs) 
leads  to  some  Important  consequences  for  the  electrical  properties  of 
the  mesosphere . 

Variations  In  the  altitude  distributions  of  the  ionization  rate  q 

or  the  loss  characteristic  P,  as  characterized  by  the  scale  factors  Z^ 

and  Z. ,  produce  a  variation  in  the  electron-density  structure  and  in  the 
P 

electrical  properties  of  the  mesosphere. 

We  examine,  first,  a  limiting  case  where  q  Is  uniform  In  altitude 

(l.e.,  Z  ■  <»),  and,  second,  the  variation  of  Z.  or  Z  In  non-uniform 
'•q'''  '  P<1 

Ionization  such  as  cosmic-ray  activity  at  night. 

Figure  1  shows  how  the  electron- loss  parameter  P  depends  on  altitude 
(due  to  such  variations  as  air  density  and  composition).  With  uniform 
ionization  (l.e.,  Z^  ■  «°),  a  doubling  of  the  altitude  scale-factor  can 
cause  a  sixfold  Increase  in  electron  density  (for  example,  at  an  altitude 
of  40  km,  between  Z^  *  5  km  and  Z^  "  10  km) . 

Figure  2  shows  another  class  of  effects  due  to  the  variation  of 
one  altitude  scale-factor  relative  to  another.  We  see  a  non-uniform 
production  of  Ionization  with  a  single  exponential  of  characteristic 
scale  Z  ■  6  km.  Ihls  value  Is  selected  as  approximately  equal  to  the 
atmospheric-density  scale  height;  It  thus  represents  Ionization  production 
that  Is  proportional  to  atmspherlc  density  and  whose  single  source 


Relotive  elsctron  density 


-7- 


Altitude  in  km,  z 

Fig.  1  —  Dependence  of  electron  density  on  voriotion 
of  electron— ‘loss  factor  (or  its  equivalent, see  text) 
for  limiting  cose  of  electron —production  foctor  Zq*a) 
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Is  ultra-hlgh'snergy  radiation,  e.g.,  cosmic  radiation  at  nlg^t.  We 
see  that  abrupt  changes  In  the  height-distribution  of  the  electron 
density  can  occur  for  changes  In  the  atmosphere  that  produce  quite  small 

variations  of  the  altitude  scale-factor  for  electron  loss,  Z. . 

P 

The  Z  factor  for  the  electron  density  has  a  symmetrical  form 
except  for  signs  In  the  Z^  and  Z^  components.  Correspondingly, 
therefore,  "trigger  conditions"  or  circumstances  for  abrupt  changes 
In  the  mesosphere's  electron  density  can  also  exist  for  small  variations 
In  Z^,  l.e.,  In  the  altitude  dependence  of  the  Ionizing  mechanism.  An 
abrupt  change  In  the  electron  density  distribution  of  the  mesosphere 
can  result,  for  exanq;>le,  from  variations  In  the  energy^  distribution  of 
the  Incoming  Ionizing  radiation.  This  effect  Is  seen  In  excitation  by 
corpuscular  radiation,  treated  by  Reid  (1961) .  He  finds  that  hlg^-lntenslty 
fluxes  of  lower-energy  solar  particles  produce  an  abrupt  transition,  both 
In  electron  production  at  about  60  km  In  the  middle  mesosphere,  and  In 
electron  density  at  about  80  km  In  the  upper  mesosphere. 
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IV .  EFFECTS  OH  THE  EUCTRICAL  MORPHOLOGY  OF  THE  MESOSPHERE 

Some  noteworthy  consequences  for  the  mesosphere's  electrical 
properties  follow  from  the  above  considerations.  They  can  be  seen 
by  examining  the  complex  index  of  refraction,  which,  for  most  of  the 
mesospheric  region  under  consideration,  can  be  approximated  by  the 
expression 

M  -  |^u)p/(vu))^J  exp  (-jTT/4)  , 

as  Indicated  in  Goldberg  (1962) • 

In  terms  of  the  electron  density,  with 

uj2  -  3180  Ug  -  , 

P 

M  then  takes  the  form 

M  -  (K/(b)^  exp  (-jTT/4)  , 

Using  the  exponential  expression  for  the  electron  density  variation 
with  altitude, 

“e  "  “eo  ®*P<'o/2qp)  exp(z/Zqp)  , 

and  an  exponential  model  of  the  collision  frequency  based  on  the  values 
of  Hlcolet  (1959), 

V  -  exp(z^/Zv)  exp(-a/Z^)  , 

we  see  that 

M  -  (K/«b)^  ^ 

From  this  expression  we  can  obtain  a  composite  altitude  scale-factor 
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given  by 

■  <%>'‘  *  > 

characterizes  the  rate  at  tAilch  the  mesospheric  Index  of  refraction 
M  varies  with  altitude,  since 

M  St  (K/o))^  C**eo^**ov^^  exp(z/Zjj)  exp(-jTT/4)  , 

where 

■  ‘ov/^v>]  . 

One  physically  significant  condition  Is  the  limiting  value  of 
Z^.  Thus  when  electron  production  varies  more  slowly  in  altitude 
than  does  electron  loss  (or  Its  composite  equivalent  In  the  expression 
above),  Z^  In 

Zp  Z^j 

must  always  be  less  than  2Z^,  idilch  Is  about  2x6.1*  12. 2  hm.  The 
electrical  dissipative  process  In  the  atmosphere  (characterized  by  v) 
not  only  provides  local  damping  to  electromagnetic  effects,  but  also 
limits  the  rate  at  which  such  effects  can  vary  with  altitude.  Because, 

In  these  regions  of  the  atmosphere,  v  depends  principally  on  the 
atsiospherlc  density,  variations  in  the  limiting  value  of  the  local 
refractive  Index  will  depend  principally  on  variations  of  the 
atsBOspherlc  density. 

Ion  production  by  varying  mechanlssis  and  under  differing  atmospheric 


conditions  can  cause  the  mesosphere's  electrical  properties  to  have 
many  variations  with  altitude,  soma  with  limiting  conditions.  Table  1 
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Table  1 

ALTITUDE  SCALE-FACTORS  FOR  DIFFERING  CONDITIONS  OF 
lONIZATIOi  IN  THE  MESOSFHERE 


Case 

Reference 

^n 

fkml 

Si 

fkml 

Daytime;  60-70  Ion 

Fejer  and  Vice 

8.2 

7.2 

Daytime;  70-80  km 

Fejer  and  Vice 

13.8 

8.2 

Nighttime;  65-80  km; 
cosmic  radiation 

Moler 

5.5 

5.6 

Daytime;  50-70  km; 
cosmic  radiation  and 
Lyman  a  radiation 

Holer 

4.2 

4.8 

Daytime;  30-60  km;  solar 
protons 

Reid 

4.1 

4.8 

Nighttime;  70-90  km; 
solar  protons 

Reid 

-10.1 

3.0 

Summer  day;  55-65  km 

Landmark  and  Lied 

5.2 

5.4 

Summer  day;  65-70  km 

1 

Landmark  and  Lied 

-3.5 

-17.1 
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■unoarices  ■oma  of  the  states  of  Ionisation  that  have  been  treated  In 
the  literature.  For  the  electron-density  factor,  Z^,  and  the  refractlvlty 
parameter,  Z^,  It  gives  values  computed  for  those  portions  of  the  n^ 
altitude  profile  that  can  be  represented  approxlisately  In  terms  of 
s  single  exponential  characteristic. 

Figure  3  shows  the  altitude  scaling-factor  for  mesosphere 
refractlvlty,  Z^.  In  this  case,  the  collision-frequency  altitude 
acallng-factor  Z^  Is  fixed  at  6  km. 

Note  that  the  mesosidiere  can  have  three  distinctive  conditions  of 

Z  .  These  are: 
n 

1.  Electron  density  Increasing  with  altitude  (upper  rl^t  In 
Fig.  3).  In  this  condition,  the  refractive  Index  Increases  with 
altitude  but  does  not  vary  more  slowly  than  Z^  ■  12  km. 

2.  Negative  variations  of  electron  density  with  altitude,  «diere 
Z^  ranges  from  sero  to  -6  km  (lower  left  In  Fig.  3).  The  refractlvlty 
rapidly  falls  off  with  Increasing  altitude  for  values  of  Z^  lying 
between  zero  and  2  km.  For  ^  3  km,  the  refractlvlty  changes  slowly 
with  altitude,  with  Z^  varying  from  -12  km  (for  Z^  ■■  -  3  km)  to  Infinity 
(for  Z^  ■  -  6  km) . 

3.  Electron  density  falling  off  with  altitude  at  a  rate  of  Z  >  -6  km 

It 

(upper  left  In  Fig.  3).  This  leads  to  gradual  upward  Increases  of 
refractlvlty. 

Variations  with  tlise  In  the  excitation  mechanism  producing  q 
(e.g.,  variations  In  flux  or  Intensity  of  cosmic  radiation,  solar  proton 
boodbardaent ,  or  other  processes  affecting  the  state  of  Ionisation)  will 
cause  corresponding  fluctuations  In  Z^.  This  will  lead  to  tlsM  variation 


Slow  upword 
incroosM  iiv 
refractivi>^ 


40  Z|,  in  km 


■50  -40  -30  -20  -10 


Z,,*  -12  km - 


roftronet  volvt 


Ropid  upward  incrtosM 
in  rtfroetivity  — i 

12  km  limit/ 


20  30  40  50 

Zn  in  km 

Rapid  dawnward  incrtatot 
in  rafractivity 

Slow  downword  incraom 
in  rtfroetivity 


Fig.  3  —  Four  conditions  of  Zm  (the  refroctivity  scale 
structure)  in  the  mesosphere  for  various  conditions 
of  variation  of  Zn  — electron  density  with  altitude 
(case  of  collision— frequency  altitude  variation 

Zj/ «  6  km) 
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o£  th«  electron  density  In  the  neeoephere  end  in  the  etructure  lectors 

Z  or  Z.  ,  resulting  in  a  "isodulation"  of  the  electrical  structure 
n  Pq 

factor  Z^.* 

The  variations  can  occur  for  any  of  the  four  conditions  of 
mesosidierlc  ionization  shown  in  Fig.  3.  However,  the  Dost  sensitive 
atnospheric  conditions  (causing  great  responsiveness  to  changes  in  the 
Dodes  of  excitation)  occur  when  the  atmosidiere  is  in  those  states 
•unaerized  in  Thble  2. 

Vltether  or  not  the  value  of  Z  varies  considerably  with  altitude 

<1 

dependa  on  the  flux  and  energy  distribution  of  the  ionizing  radiations. 
The  nesos^ere's  electrical  morphology  can  then  be  described  in  terms 
of  segments  of  exponential  variations  of  the  parameter  M.  Under  complex 
conditions,  the  varied  electrical  structuring  depicted  in  the 
diagram  can  occur  sioniltaneously  at  various  altitudes.  This  could  happen 
uhen  the  mesosphere  is  excited  by  highly  enhanced  ultraviolet  or 
x-radiation  or  by  solar  particles  during  solar  stonu.  The  diagram 
Indicates  that,  at  that  time,  shelves  or  strata  could  be  a  prominent 
BMSospherlc  structure. 
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Table  2 

COHDITIONS  OF  MESOSPHERE  FOR  SENSITIVE  RESPONSE 
IN  ELECTRICAL  STRUCTURE  M  BY  VARUTIOM 
OF  lONlZATlO)  DISTRIBUTION  WITH  ALTITUDE 


Condition 

Altitude  Variation 
of  Electron  Density 

Mesosphere  Response 

1 

0  <  <  6  km 

M  Increases  rapidly  with  altitude 
for  small  increases  of  Z  and  is 
positive 

2 

-2  <  <  0  km 

M  decreases  rapidly  with  altitude 
for  small  decreases  In  Z  and  is 
negative 

3 

-6  <  Z  <  -3  km 
a 

M  becomes  'diffuse  downwards'  (i.e., 

Z^  is  large  and  negative) »  increasing 
slowly  with  decreasing  altitude;  large 
positive  gradient  changes  in  M  occur 
for  small  decreases  in  Z 

n 

4 

-15  <  Z  <  -6  km 
n 

M  becomes  'diffuse  upwards'  (i.e., 

Z^  is  large  and  positive),  ircreaslng 

slowly  with  increasing  altitude; 
large  negative  gradient  changes  in 

M  occur  for  small  decreases  in  Z 

n 
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